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Recent experiments on quantum dots in the Coulomb Blockade regime have shown how adding 
successive electrons into a dot modifies the energy spectrum and the wave functions of the 
electrons already present in the dot. Using a microscopic model, we study the importance 
of electron-electron interaction on these "scrambling" effects. We compute the Hartree-Fock 
single particle properties as function of the number p of added electrons. We define parametric 
correlation functions that characterize the scrambling properties of the Hartree-Fock wave 
functions and energy spectra. We find that each of these correlation functions exhibit a 
universal behavior in terms of the ratio p/p* where p* is a characteristic number that decreases 
with increasing either the interaction strength, the disorder strength or the system size. 



1 Introduction 

At very low temperature, by applying a gate voltage Vg to a, Quantum Dot (QD) weakly con- 
nected to leads by tunnel contacts, one obtains a series of Coulomb blockade conductance 
peaks of height Gn at positions . The peak positions are directly related to differ- 
ence En+i — E]\f between the ground state energies of the dot with and A -|- 1 electrons (in 
zero bias). The heights Gn are related to the overlap of the | A > and | A -|- 1 > eround state 
wave functions at the contacts 0. During the last few years, many experimentsHS have shown 
that the statistical properties of the successive peak positions could not be described within 
the framework of the random matrix theory (RMT). For the addition spectrum (e.g. statistical 
properties of the differences A2 = Kj^^^ ~^g^)i qualitative agreement with experimental results 
can only be obtained by taking into account fluctuations of the charging energy that originate 
from the spatial fluctuations of wave function amplitudes i'Bi. Although such fluctuations of 
wave functions already exist in RMTq, the resulting charging energy fluctuations are paramet- 
rically to small if one ignores the reorganization or "scrambling" iof the wave functions and of 
the spectrum) induced by the addition of the successive electrons BS. In fact, in a recent exper- 
imental study of the peak to peak heights correlation, Patel et al. i have unambiguously shown 
that this reorganization is important and could essentially be attributed to electron-electron 
interaction effects. Recently, Alhassid et al. @ have confirmed the importance of scrambling us- 
ing a phenomenological model based on Gaussian Random Matrix Process (GRMP^, and have 
obtained good agreement with Patel et al. experiment. 

In the present work based on a microscopic model of a quantum dot, we present a de- 
tailed numerical study that shows directly the importance of the electron-electron interaction 
on the scrambling of the single particle properties obtained within a self consistent Hartree-Fock 
(SCHF) method. 

2 Microscopic model 

We model a two dimensional disordered dot of size A = x Ly hy the following Hamiltonian: 

U ^ 

H = J2 Wicjci - t ^ [c|cj + h.c] + ^ij^i^V (1) 

i <ij> ')i=l 

where < ij > denotes nearest neighbors on a square lattice of size A = Lx x Ly and t = 1. cl 
and Ci are creation and annihilation operators for an electron (spinless) on site i and Ui = c[ci. 



The random on-site disorder potential Wi are uniformly distributed over [— ^,^]. Because 
of the disorder, the electronic charge density exhibits fluctuations in space. We consider that 
the natural bare interaction between these charge fluctuations has a long range Coulomb form 

^ii = uT^n- We choose a torus geometry to prevent from electrostatic effects that tend to 

I * j\ 

accumulate the charge on the surface. The distance |rj — rj\ is calculated from the shortest path 
relating the two sites i and j on the torus. 



3 Hartree Fock picture of scrambling 



We treat interaction effects in the self-consistent Hartree-Fock approximation. We thus assume 
that the ground state I^'q' > of a dot with particles is a Slater determinant of N effective 
single particle states i/^^ (r^) with associated quasi-energy e^. These states Tp^ (rj) are in fact 
the self-consistent eigenstates of an effective one-particle Hamiltonian = J^i'^i^'^ ^Ici — 
J^Kijyitij^ (^l^^j + h.c.]. The effective onsite potential w^^^ = Wi + UY^jMij < ^Q^|nj|^'Q^ > is 
the sum of the disorder potential and the Hartree contribution. Similarly the effective hopping 
^i/^ = i(J<ij> + UMij < ^'Q^|c|cj|^o^ > includes the Fock exchange term. Self consistency is 
obtained with an iterative procedure. 

The self consistent states (rj) are linear combinations of the non interacting 'ip^{ri). 
However the ij^^ (rj) depend explicitly on interaction strength and particle number and not 
only on disorder strength and system size as the Vo(ri). Therefore, by adding one by one p 
electrons into a dot with initially N particles, the effective one body Hamiltonian is changed 
into another Hamiltonian h^^'P with the corresponding il>^^P{Yi) and e^^^ • the Hartree- 
Fock approximation, one can always write h^~^P = + y^'^+P , such that the scrambling 
effects of the p additional electrons is described by an effective one body term y^'^+P, that 
is neither diagonal in the position basis nor in the V',^ basis of the dot with particles. In 
fact this possibility of writing /i^+P = + y-^'^+P allows us to make a parallel with the 
phenomenological model of Alhassid and MalhotraB. In their work, the N+p interacting particles 
system is described by an effective one body Hamiltonian /i^+P = cos {p5x) + hi sin {pSx) 
in which Hq , hi are two random Gaussian matrices and 5x is the parameter that control the 
scrambling strength. Their effective scrambling potential y^'^+P is thus given by y^'^+P ~ 
/io(cos {p5x) — 1) + hi sm{p5x). 



4 Numerical results 



4-1 Scrambling properties of SCHF single particle wave functions 

In order to characterize the scrambling of the SCHF single particle wave functions, we have 
computed the following parametric correlation function: 



dr +^(r)V^*(r) 
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(2) 



in which t-t-t means averaging on disorder configurations and on the A^ lowest single particle 
energy levels a. The scrambling of wavefunctions is characterized by a decrease of C'^(p) with 
increasing p. A similar quantity has already been studied in the context of GRMPB. We have 
studied in detail the dependences of C^{p) on the disorder strength W (or conductance 5), 
the interaction strength U and the system size A. The ranges of parameters considered are 
ly = 4 ^ 10; [/ = 0.5 ^ 3, A = 6 X 6 ^ 11 X 10 sites. The initial values of A^ are chosen so that 
the band is quarter filled v = N/A ~ 1/4. For each set of parameters we have averaged over 500 
distinct disorder configurations. The results obtained for C^{p) are summarized on Figure la. 
This figure shows a quite interesting result, that the correlation function C^{p) can be recast 



into a one parameter scaling form C^{p) = f^{p/pip) where depends on the three parameters 
W,U,A. In the range of parameters studied we further obtain the fohowing approximate fitting 
form: 

l/p^ =cU^A^W^[l + dAU\ 

where c ~ 8.10"^ and cLa = 0.03; 0.05; 0.1; 0.15; 0.2 for A = 6 x 6; 6 x 7; 7 x 8; 8 x 9; 9 x 10 
respectively. For a given p, C^{p) decreases when either W,U ox A is increased. Thus the larger 
the system size, the stronger is the scrambling effect of a single additional electron. Let us 
compare the scaling (3) with the scaling obtained by Alhassid and Malhotra in the context of 
the GRMP mo dell. The behavior found in this model for the correlation function is C^^*^^(p) ~ 
[1 + (^)^]~^ with 1/p^ — 6x^/M for gaussian orthogonal matrices of size M x M. The difference 

between the functional form C^{p) and C^^'^^^ (pgrmp) seems to prevent any simple mapping 
between them. Nevertheless, if we consider only the tail of the two forms, we can identify the 
microscopic value p with Pq^mp ■ then obtain the following relation between the microscopic 
parameters and the parameters of the GRMP model : ((5x\/M)^^/^ ~ c U'^A^W^{1 + dA U"^)- 




Figure 1: (a) C^{p) as a function of the scaling variable p/p^. (b) Ct{p) as a function of the scaling variable 
p/Pf For a given disorder W and system area A, we plot all the data points corresponding to the entire range 

p — {1, 2, .., 10} and interaction range U = {0.5, 1, .., 3}. 

4-2 Scrambling properties of SCHF single particle energy spectra 

In order to characterize the scrambling of the SCHF single particle energy spectrum we have 
computed the following parametric correlation function: 

N+p 

where = e^+i ~ is the spacing between two consecutive levels {s^ is normalized to the 
mean level spacing at this energy) . Since the one particle density of state has a deep around the 
(average) Fermi level, we have used all a corresponding to energy levels (empty or occupied) 
neither in the this deep nor at the band edges. As it is illustrated on Figure lb, we have found 



that Ce{p) can also be recast into a one parameter scaling form C^ip) = feip/Pe) where 
depends on the parameters W, U, A. Quantitatively we obtain the following approximate fitting 
form: 

where c ~ 5.10~^. The quantitative behavior of p^, in terms of W, A, U, is different from that of 
p^. This contrasts with GRMP-like models in which necessarily p^ = Pip- For the microscopic 
model, we believe that the difference found between p^ and p^ reflects the distinct statistical 
behavior of the diagonal and the off-diagonal elements of y^'^+P in the ip^ basis (within 
a perturbative approach, the SCHF single particle energy level position depends on both the 
diagonal and off-diagonal terms. In contrast the reorganization of SCHF wave functions depends 
only on the off-diagonal elements). 

5 Conclusion 

Based on a microscopic model of a disordered quantum dot, we have presented detailed numerical 
results that show the importance of electron-electron interaction on the scrambling of the single 
particle properties obtained within a self consistent Hartree-Fock method. As main result of this 
study, we have found one-parameter scaling properties for the two quantities (p) , (p) which 
describe respectively the reorganization of the one particle wave functions and of the energy 
spectrum when p particles are added to the system. The quantitative analysis of this one pa- 
rameter scaling properties cannot be described by effective GRMP models considered until now 
01. This suggests that new specific random matrix models are necessary to describe microscopic 
disordered models of Coulomb interacting particles within the Hartree-Fock approximation. 

References 

1. L. P. Kouwenhoven et al., in Mesoscopic Electron Transport, Edited by L.L. Sohn, L.P. 
Kouwenhoven, and G. Schoen (NATO Series, Kluwer, Dordrecht, 1997) 

2. U. Sivan et al, Phys. Rev. Lett. 77, 1123 (1996), 

3. F. Simmel, T. Heinzel and D. A. Wharam, EuroPhys. Lett. 38, 123 (1997), S. R. Patel 
et al, Phys. Rev. Lett. 80, 4522 (1998), F. Simmel et al, Phys. Rev. B 59, 10441 (1999) 

4. Ya. M. Blanter, A. D. Mirhn and B. A. Muzykantskii, Phys. Rev. Lett. 78, 2449 (1997), 
R. O. Vallejos, C. H. Lewenkopf and E. R. Mucciolo, Phys. Rev. Lett. 81, 677 (1998) 

5. R. Berkovits, Phys. Rev. Lett. 81, 2128 (1998), S. Levit and D. Orgad, Phys. Rev. B 60, 
5449 (1999), P. N. Walker, G. Montambaux and Y. Gefen, Phys. Rev. B 60, 2541 (1999), 
A. Cohen, K. Richter and R. Berkovits, Phys. Rev. B 60, 2535 (1999), L. Bond and R. 
Berkovits, Europhys. Lett. 47, 708 (1999), K-H Ahn, K. Richter and LH Lee, Phys. Rev. 
Lett. 83, 4144 (2000) 

6. S. R. Patel et al, Phys. Rev. Lett. 81, 5900 (1998), 

7. Y. Alhassid and S. Malhotra, Phys. Rev. B 60, R16315 (1999) 

8. Y. Alhassid and H. Attias, Phys. Rev. Lett. 74, 4635 (1995), ibidFhys. Rev. E 52, 4776 
(1995) 



